Berry phase in graphene: a semiclassical perspective 
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We derive a semiclassical expression for the Green's function in graphene, in which the presence 
of a semiclassical phase is made apparent. The relationship between this semiclassical phase and 
the adiabatic Berry phase, usually referred to in this context, is discussed. These phases coincide 
for the perfectly linear Dirac dispersion relation. They differ however when a gap is opened at 
the Dirac point. We furthermore present several applications of our semiclassical formalism. In 
particular we provide, for various configurations, a semiclassical derivation of the electron's Landau 
levels, illustrating the role of the semiclassical "Berry-like" phase. 
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I. INTRODUCTION 

Graphene [l!, , a two-dimensional carbon based ma- 
terial forming a honeycomb lattice, has attracted a lot of 
attention since its experimental isolation has been proved 
possible d, It is a gapless semiconductor in which, 
near half filling, electrons behave like massless Dirac par- 
ticles, obeying a linear dispersion relation. Among the 
unusual properties of this two-dimensional carbon mate- 
rial stand out very distinctive quantum Hall properties, 
and in particular the -^n dependence of the energy in 
terms of the Landau level number n, and the existence 
of a Landau level with zero energy, which is associated 
with the presence of a Berry phase [l|, Q . 

The existence of this Berry phase and its implications 
for the Landau levels have been discussed in many places 
in different contexts (see e.g. [1,H,0|)- The direct connec- 
tion between the Berry phase and the observable quanti- 
ties under discussion is however not always as transparent 
as one may wish, and situations where, either because of 
disorder, or because one would like to confine the elec- 
trons into a finite region of space, a position dependent 
electrostatic potential or mass term is introduced, are 
usually not addressed. 

The aim in this paper is to revisit this question of 
Berry phase in graphene within a semiclassical, and more 
specifically semiclassical Green's function, perspective. 
For sake of clarity, our emphasis in this present work 
will be more in providing this new point of view, and we 
shall therefore mainly illustrate it with the discussion of 
the standard problem of the Landau levels of electrons 
in a perpendicular and uniform magnetic field. Even in 
this familiar framework, we shall see however that our 
semiclassical approach makes it possible to address some 
non-trivial questions, such as the role of the Berry phase 
in situations for which a small mass term has to be in- 
cluded, opening in this way a gap at the Dirac point. 

This article is therefore organized as follows. In sec- 
tionini we derive, following closely the formalism of Boltc 
and Keppeler [Sj, the expression for the semiclassical 
Green's function in graphene. In particular we discuss 
in details the origin of the term corresponding to the 



Berry phase. These results are extended in section Hill to 
a bilayer of graphene. We furthermore provide both for 
the monolayer and the bilayer cases the expression of the 
Gutzwiller trace formula for the semiclassical density of 
states, valid when classical periodic orbits are isolated. 
As an illustration of the Green's function formalism, we 
then apply it in section HVl to the computation of Landau 
levels for a graphene sheet in constant magnetic field. We 
will see in particular that the modifications brought in by, 
for instance, trigonal warping, are easily included within 
our semiclassical formalism. We then come back in sec- 
tion |V] to the discussion of the relationship between the 
semiclassical "Berry-like" phase obtained in our approach 
and the adiabatic Berry phase [9,] usually discussed in this 
context. 



II. SEMICLASSICAL GREEN'S FUNCTION 
FOR GRAPHENE 

Starting from a tight-binding nearest neighbor model, 
the graphene Hamiltonian at low energies can be ob- 
tained by expanding the momentum near the Dirac 
points K and K' of the Brillouin zone. For pure 
graphene, one obtains in this way in momentum repre- 
sentation (10,. .11.. .12... 13.] 



<yyPy) 



apx - ipy 
apx + ipy 



(1) 

where the matrix structure originates from the exis- 
tence of two sub-lattices (denoted A and B below) in 
the graphene honeycomb structure. In this equation, 
Vp = Sta/(2h) is the Fermi velocity, with t the hopping 
parameter and a the lattice constant, a is the valley index 
{a = ±1) labelling the two inequivalent points K and K' 
in the Brillouin zone (not to be confused with the sub- 
lattice index), p is the momentum measured from these 
points, and cr^j.y are Pauli matrices. This linear approxi- 
mation to the graphene Hamiltonian will be valid as long 
as the condition |p| H/a is fulfilled. 

We are interested here in a more general situation than 
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the one of pure graphene, and would like to consider the 
case where, because of either disorder or the need to con- 
fine the electrons in some part of the graphene sheet, an 
electrostatic potential U{r) and/or a [possibly position 
dependent] mass m(r) have to be taken into account. 
We will not consider however tunneling contributions re- 
lated to the Klein paradox, or boundary effects that may 
occur at the (zigzag, armchair, or generic) edges of the 
graphene sample. The graphene Hamiltonian then takes 
the more general form 

Hg = v^iaa^K^ + aylly) + C/(r).l2 + m{r)vla, , (2) 

in which the magnetic field B(r) = V x A(r) (if any) is 
taken into account by the Peierls substitution 

p ^ n = p + eA(r) , (3) 

with A(r) the vector potential and p = —ih-^. 

For this problem, the Green's function G(r",r') is ac- 
tually a 2 X 2 matrix defined by the differential equation 



{E.U~Ug)G{v",v'-E)^5{v" -v').U 



(4) 



(where Tig is applied to the variable r"). To obtain a 
semiclassical solution of this equation, we shall proceed 
in two steps. First, assuming r" is far from the source 
location r', we solve semiclassically (i.e. in the WKB 
approximation) the Schrodinger equation 



{E.l2~Hg)G^Q 



(5) 



In a second stage we match this general solution to the 
exact Green's function of the "free" (i.e. with constant 
potential and mass) problem, valid near the singularity 
r'. We proceed now with this derivation. 



A. Far from the singularity: the WKB 
approximation 

Following fsl] , we seek a semiclassical solution of eq. ^ 
with G of the form 



G(r",r';£;) =r(r",r')exp 



-5(r",r') 



(6) 



where F is a 2x2 matrix. To lighten the notation, we drop 
for now the explicit dependence in the source position r'. 
Inserting ^ into ^ and expanding in h the resulting 
expression, we obtain at order 0{hP) 



i?.l2-ff(^,r"))r(r")=0, (7) 



and at order 0{h^) 

f .Ar(.",.„,(„.^ + <,,|,)r(r")^o.(8) 

where H{p,r) is the classical symbol associated with the 
quantum Hamiltonian Tig. 



This classical Hamiltonian can be diagonalized, with 
the eigenvalues 



H^{p,r) = U{r) ± y/m^{r)vj+vim 



(9) 



and the corresponding normalized eigenvectors (p, r) 
(whose explicit expressions are given in appendix 
Writing the matrix r(r") as [F±(^, r") • f ±(r")], with 
a 1 X 2 matrix, the order hP equation becomes 



(10) 



where the ± sign must be taken according to the sign of 
E~U{r"). 

Eq. pO|) is the usual scalar Hamilton-Jacobi equation, 
which can be solved by the method of characteristics 
p^ . This amounts to constructing a 2-dimensional La- 
grangian manifold C (in the 3-dimensional energy surface 
in phase space) built as a 1-parameter family of trajec- 
tories following the classical equations of motion 



dH 



± 



dp 

P = ^ P,r) ■ 

or 

Given any such manifold, the action S{r") = pdr, 
where the integral is taken on an arbitrary path on £, is 
a solution of pO)) . 

The specific Lagrangian manifold that will correspond 
to the proper boundary conditions for G(r", r') near the 
source r' is the one obtained from the trajectories leaving 
r' with an arbitrary initial momentum p' at energy E: 



£± ={(p(t),r(i)),t e [0,oo), 

such that r(0) r', (p(0), r(0)) = E} 



(11) 



(each point on the manifold is therefore parameterized 
by the time t and the initial momentum p(0)). The cor- 
responding action can then be expressed as 



5±(r",r') p-rdt 



(12) 



along a trajectory (p(i),r(i)) joining r' to r" at energy 
E. 

Having obtained a solution of the 0{hP) equation, the 
prefactor F is then determined by the 0{h^) equation ([H]), 
which, after multiplication on the left by ^^^(^77-, r"), 
can be expressed as DF* = 0, where 

□ =ry±t(^ r")^ AV±(^ r") 

V dv" ' ' dp dv") ^dr" ' ' ■ 



The operator □ can be decomposed as □ = □( 
with 



1) 



□ 



(2) 



dp 



d_ 



(13) 
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and 



□ 



(2) 



9p \ dv" 



(14) 



Noting that first order perturbation theory impfies 
V^'^[dH/dY>)V^ = (9i/±/ap), one has straightfor- 
wardly that 



and that 



(15) 



for any system of coordinates) and /i^ is the Masfov in- 
dex counting the (algebraic) number of caustic points. 
Writing 
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we obtain that 



and therefore S*(<) — exp {i£,sc) S*(t = 0), with 



op 



(Note here that with respect to spatial derivation, = 
ff±(r") = H^{{dS^/dr")y)). One recovers in this 
way, for the real part of □, the usual expression valid for 
a scalar quantum system [ij], which is expected since it 
basically expresses the conservation of probability. 

The imaginary part of (0(2)) however is not con- 
strained by such a conservation law, as it affects only the 
phase of F, but encodes information about the adiabatic 
variation of the eigenvector along the followed trajec- 
tory. It needs therefore to be computed from the explicit 
expressions of the eigenvector and eigenvalues of (p, r). 
The details of the algebra are given in appendix [X] One 
obtains 



\ 9p dr" 



1 d dH^ 
2d?'' dp 



f± = 



(17) 



with 



2{E-U{v")) 

_d 
dr 



eB 



nxj^{m{r")vl~U{r")) 



(18) 



m{r")vl+E-U{r") 



is the unit vector in the direction perpendicular to 
the graphene sheet). 

In the absence of the complex term iM, the scalar 



transport equation y 
the usual solution 



1 d dH^ 



2 dr" ■ dp 



r 



± _ 



c 



J^(r",r') = 



J±(r",r')| 



has 



(19) 



(20) 



where ry and r± are the coordinates parallel and trans- 
verse to the trajectory (actually Eq. (|20p remains valid 



M^{p{t'),r{t'))dt' 



(21) 



Summing the contributions corresponding to different or- 
bits j joining r' to r" we get 

G(r",r';i?)= J2 7^ V-/(r")S,±(t = 0) 

j:r'^r" 

exp f^5f(r",r') - ^ Mf{p{t') Mt'))dA , 



(22) 



where T^/(r") = V^{dSf /dr" ,t") (and therefore de- 
pends not only on r" but also on the final momentum 
p" of the trajectory j). 

The semiclassical phase ^sc Eq. (PT|) is the analog, in 
our context, of a Berry phase |9j]. In the same way, it has 
its origin in the adiabatic change of the eigenvectors of 
the "internal degree of freedom" Hamiltonian if (p(r), r) 
along the classical paths contributing to the semiclassical 
Green's function. Furthermore, in some circumstances, 
^sc exactly corresponds to the genuine Berry phase ^ad de- 
fined for the adiabatic motion along the trajectory. This 
will be the case in particular for "pure" (i.e. without 
mass ter m) g raphene. In general, however, £^sc and ^ad 
differ [l^, [l^. We will come back to this point in sec- 
tion |Vl and in particular clarify the question of which of 
the two phases is relevant for the Landau levels. 



B. Matching to the exact solution near the source 

Sufficiently close, on the classical scale, to the source 
r', we can neglect the variation of the various potentials 
and of the mass, i.e. assume U{v) = Uq, m(r) = mo and 
A(r) = 0. In this case we have the expression for the 
exact retarded Green's function: 



Q — ( Gaa Gab 
\Gba Gbb 



with 
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GAAir",r',E + ie) - 
GABir",r',E + ie) = 



* — — " 



-ioLcj) 



C 



(23) 
(24) 



r 



Gaa(™o -"to), G 



__BA 



is the phase of 



and Gbb - 

Here C = ^/{E + ie - Uo) 
Px + and Hq and i?i are Hankel functions of order 
and 1. Asymptotically, as |r" — r'| +oo, Gaa and 
Gab take the form 



Gaa 



.-^^-l+E^U, f2em^"~r,^,) (25) 



2 e 



4(;ifp)2 



i(fc|r"-r'|-f ) 



(26) 



with = -L^^E-Uo)^-m'o<^ IpI- 

Let us assume E — Uq > 0, so that semiclassically 
we consider the positive eigenspace . We note first 
that, in the free case considered here, the choice of the 
Lagrangian manifold given by pT|) corresponds to the 
action S'+(r", r') = |p| • |r" - r'| and to 

so that, as anticipated, the expression ([22]) matches the 
asymptotic expressions ([25l) -(|26 l) . provided one chooses 



G = 



and 



E+(t = 0)^F+t(^,r'). 



(27) 



The asymptotic expressions (I^St - (|26l) are valid as soon 
as |r" — r'l is larger than a few Fermi wavelengths, which 
can still correspond to a distance short on the classical 
scale, and therefore such that the free Green's function is 
a good approximation. We can therefore use this match- 
ing condition to fix the prefactors I]+(t = 0) and C in 
the generic case, obtaining finally 

G(r",r';i?)= ^ ^ 



exp 



3 



III. 



V^{v")-V^\v'). 

(28) 



BILAYER GRAPHENE AND GUTZWILLER 
TRACE FORMULAE 



We turn now to a few extensions of the result derived 
in section |TT1 We start with a generalization to the bi- 



layer graphene case, and then briefly discuss the resulting 
Gutzwiller trace formulae for the density of states, valid 
when classical periodic orbits are isolated in phase space 
(i.e., generically, for chaotic systems). 



A. Semiclassical Green's function for the bilayer 

case 

The bilayer graphene Hamiltonian can be written at 
low energy as [l7| 



-T/U _ 



1 



2ra* 





{Px + iPy 



{Px - iPyf 
|2 



(29) 



with m* = 7i/(2u^), where 71 is the intra-layer coupling 
parameter. As before, we would like to include electric 
or magnetic fields, as well as a possibly position depen- 
dent mass term. We therefore consider the more general 
Hamiltonian 



Hu = C/(r).l2 + ni{v)vla, + H^,(p 



n) 



(30) 



Following the same approach as above, one obtains the 
semiclassical Green's function in the form Eqs. 
except for a different expression of the classical Hamilto- 
nian eigenenergies 



U{y) ± Wm(r)2i;4 + 



/ IP 



\2m* 



and of the semiclassical ( "Berry-like" ) phase term 




m(r")^u^ 



±eB - 



[E - U{y")Y 

1 n X d[m{r")vl - C/(r")]/ar" 

2 m{v")vl + E- U{r") 



(31) 



In the free case (m(r) = mo, U{r) = Ua), the exact 
Green's function can be shown to behave asymptotically 
as I r" — r'l +00 as 



Gaa 



'^BB ^ GAAi'mo 



movj 



+ E-U0 /2"e^(''^l''"-'"'l-T) 



G 



AB 



G 



BA 



-Too) 

2"gi(fe|r"-r'|-f ) 
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with (/) the phase of Px + ipy Matching the exact so- 
lution near the source to the semiclassical expression far 
from the source eventually gives the semiclassical Green's 
function as a sum over all trajectories j joining r' to r" 
under the classical Hamiltonian or (depending 
on the sign of (£; - t/(r'))) 



G(r",r';i?) = 



exp 



i J*' Mfdt' - i 



with given by Eq. pi)) . 



B. Trace formulae for isolated orbits 



(32) 



In particular, in Eqs. 



(ESI) 

and 



or (1321), V^{y") = V^{y') 
Pj = p'j (remember that 
Vj*(r) = Vj^{pj{r),r), so the second condition is nec- 
essary here). Therefore Tr[l/=^(r")-l/=^^(r')]|r"=r'=r = 1- 
Once this point is recognized, the calculation of posc from 
the semiclassical Green's functions is, up to the inclusion 
of the semiclassical "Berry-like" phase term <f. M^{t)dt, 
essentially the same as in the scalar case Jll [H (see 
also the particularly clear discussion in [2(|). We thus 
just quote the final results: p{E) = p'^{E) + p^{E); 
P^{E)=p^^^,{E)+pf,^{E), with 



"""p-o. ^|det(Mpo-l)| 

/ Q± _ rT^. 



(35) 



M^dt' 



One important application of the semiclassical expres- 
sions for the Green's functions is that, by taking their 
trace, one obtains a semiclassical approximation for the 
density of states p{E) = J^i^E — Ei). We have in 
mind here a quantum dot defined in a finite region of 
a graphene sheet (with the confinement imposed for in- 
stance through the mass term), and the Ei are the cor- 
responding discrete energies of the confined system. We 
will furthermore assume in this subsection the classical 
motion within the dot fully chaotic, so that all trajecto- 
ries are isolated. 

Starting from Eqs. ([28]) or ((32l) . the semiclassical den- 
sity of states can be obtained as the trace 



P{E) 



-Im j <irTr[G(r, r; £;)] 



(33) 



(where Tr is the trace on the internal structure of the 
Green's function). The smooth (Weyl) part of the density 
of states, which is associated with "zero length" orbits, 
has the usual expression pwcy\[E) = Pwoyi(^) + /^Wcyi(-^) 
with 



PWoyl(-^) 



dpdr 



5{E-H^{p,v)) 



(27rft)2 

When potential and mass terms are constant this gives 



PWeyl(^^) 



(34) 



with A the area of the graphene sheet and O the Heavi- 
side step function. 

The oscillating part posc{E) of the density of states 
can then be obtained inserting the semiclassical expres- 
sion for the Green's function in Eq. (j33p . Performing the 
integral on r in the stationary phase approximation im- 
poses that, in the semiclassical sums Eqs. (l28l) or ([32]). 
only the trajectories with identical initial and final mo- 
mentum should be kept. As a consequence, the sum 
over the index j becomes a sum over periodic orbits. 



Here M 



is the monodromy matrix, 



P- 



is the topologically invariant Maslov index [v = Q ot: 
1, depending on the sign of d^Sj/dr\, see the discussion 
in [20|), and T'ppo is the period of the primitive orbit 
(Tpo = nTppo if the orbit consists of n repetitions of the 
same path). 



IV. GRAPHENE IN A CONSTANT MAGNETIC 
FIELD 



As an illustration of the semiclassical Green's func- 
tion formalism, we consider in this section the simple 
(but useful) case of a graphene sheet immersed in a con- 
stant magnetic field, and show how some standard (and 
less standard) expressions can be easily re-obtained in 
this way. We start with the Landau levels in the mono- 
layer and the bilayer, without potential or mass term 
{U{v) — m(r) ~ 0), and assuming the low-energy ap- 
proximations Eqs. (fT|)- (P5)) of the Hamiltonian apply. 
We then study the influence of higher order corrections 
(e.g. trigonal warping) to this low-energy Hamiltonian. 
We finally consider the case where a finite mass term 
TO(r) = mo — const., is introduced. This last exam- 
ple will be used to introduce the discussion on the dis- 
tinction between the semiclassical and adiabatic Berry 
phases, with which we shall end this paper in the next 
section. 



A. Landau levels in monolayer graphene 

In the absence of confining potential or mass term, and 
with a constant magnetic field, the classical equations of 
motion in graphene are integrable and lead to cyclotronic 
motion, i.e. circular periodic orbits with period T and 



6 



radius R given in the monolayer case by 



T 



27r E 

vl eB 



R = —T 
27r 



(36) 
(37) 



Since the periodic orbits are not isolated, we cannot use 
the Gutzwiller trace formula derived in the previous sec- 
tion and we have to obtain the density of states directly 
from inserting the semiclassical expression Eq. (|28p in 
Eq. ([33| . Here however the classical dynamics is ex- 
tremely simple: there is only one primitive orbit, and 
the sum over j is actually a sum over the number of 
repetitions of this primitive circular orbit. We there- 
fore have = Etj/2, with tj — jT. Two caustics 
are furthermore traversed for each iteration of the or- 
bit, one midway through the circle, the other when the 
orbit comes back to its starting point, and the Maslov 
index is thus = 2j (note that, as discussed below, 
the last caustic should be included). Finally, the semi- 
classical "Berry-like" phase term Eq. ([TSl) reduces here 
to Mf{r{t)) = avleB/{2E) ^ const, so that 



Mf{t)dt = ajiT 



(38) 



The only technical point in this calculation is therefore 
that since, whatever the initial momentum, all trajecto- 
ries initiated in r' = r eventually return there, the final 
point r" = r is a caustic {dr'^/dp'j^ = 0) and the prefac- 
tor 1 / \/\Jj I diverges. As discussed in the appendix E of 
[2l| . this divergence can be cured using a mixed repre- 
sentation of the Green's function, i.e. by expressing the 
Green's function G(r",r') in terms of its Fourier trans- 
form G(p", y"; x' ,y') as 



Gix",y";x',y') 
1 



(39) 



A semiclassical expression for G can be derived in exactly 
the same way as for G, and leads to the same expression 



except for the transformations Sj 



J, 



-y"y'{ 



Sj — Sj 



p'^x" and 



dx" dx 



^ J, = -yV(afp^)-^ Thus 



dv" 



(40) 



which is not diverging since for the cyclotron motion 
dp"j.ldfJ^ = 1. The integral over in ([39|) becomes then 
straightforward (noting that dp'^/y" = d9, with 9 the an- 
gle made by the initial velocity with the x axis, this in- 
tegral basically provides a factor J^^ dO — 2it). Further- 
more the integration over position in Eq. p3p amounts 
to a multiplication by the area A of the graphene sheet. 



and as in section |TTTb1 Tr[V.*(r") 



±t 



(r')]|r"=r'=r = 1 



since the final and initial momenta are identical, 
therefore obtains 



\E) = 



\E\A 



+ C30 

)2 $Z 



cos 27rj 



E' 



2heBvl 



One 



(41) 



The total density of states is then p{E) — pwcy\{E) + 
p°^'^[E) with pwcyi{E) the smooth density of states 
(which is identical to the one without magnetic field) 
given by Eq. ((34l) . Using the Poisson formula, we there- 
fore have 



A 



+ 00 



with Ib = \/h/{eB) and 



En = sign [n)v p \2nheB\ 



(42) 



(43) 



We recover in this way the expression of the Landau lev- 
els as obtained in a fully quantal derivation [T^ . This 
approach furthermore provides a direct link between the 
phase §Q Mj{t)dt = ajir and the existence of a zero en- 
ergy level, as it cancels out the phase associated with the 
Maslov indices (another example of such a cancellation 
can be found in |i22i]). An alternative semiclassical deriva- 
tion of the graphene Landau levels can be obtained start- 
ing from the Dirac oscillator [23] , in the limit of massless 
carriers, provided the frequency of the oscillator is taken 
to be the cyclotronic one. 



B. Landau levels in bilayer graphene 

Considering now the bilayer case, we can proceed in 
exactly the same way as above except for two differences. 
First, the period T and radius R are now given by 



2tt m* 

1 = = ZTT — — ■ 

uj eB 



R = 



\E\ 



2iT'^m* 



-T . 



(44) 
(45) 



Second, the semiclassical "Berry-like" phase term 
Eq. ([nil) now reduces to Mf{r{t)) = ±eB/m* = const., 
so that 



Mj{t)dt = 2j7r 



(46) 



The Berry-like phase does not in this case compensate 
the phase associated with the Maslov index. Noting 
furthermore that, for the bilayer graphene, pwcy\{E) = 
m*A/{2Trh'^), we obtain 



A 



+00 



(47) 



B 77 — - 
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where 



(48) 



is the semiclassical approximation to the exact quantum 
values of the Landau levels, E'^^"'^"* = hw\Jn(n — 1) = 
hhjiji — 5) + 0(i). The semiclassical calculation fails 
here to account for the O(^) term. The n = and n = 1 
Landau levels, which both have zero energy, are there- 
fore not correctly described within this semiclassical ap- 
proach. However, for n > 2, the agreement between the 
semiclassical approximation and the exact result is quan- 
titatively very good. 



C. Influence of higher order corrections (in the 
parameter (a|p|/ft)) 

The next example to which we shall apply our semi- 
classical formalism is the shift of the Landau levels asso- 
ciated with deviations, for large momenta, to the linear 
a^roximation of the graphene dispersion relation Eq. ^ 

Starting from a tight-binding description of the 
graphene monolayer in which the effect of the next-to- 
nearest neighbor hopping is taken into account via the 
parameter f <C t, and expanding the resulting dispersion 
relation near the K and K' points up to third order in 
(a|p|/?i) (the reason for expanding up to third order will 
become clear below) , the resulting Hamiltonian reads (in 
the absence of electric or magnetic fields) [131 



' ' U(p) h'ip)^ 



(49) 



with H° given by Eq. IJ) and 



h\p) = _3t' + 6yi;,|p|(^|p|-2a(^)Vcos30p 

Hp) = -Vf (^("Pa; - iPyf + 2(^)2p2(ap^ + ipy) 



Keeping only terms no greater than third order in 
momentum, the eigenvalues of the associated classical 
Hamiltonian can be expressed as 



= h'ip) ± i;,|p| (1 - a^lpl COS30, 



6i 

i(^)V(3 + cos^30p)j , 



with (/)p = arctan(pj,/pa;). The anisotropic terms, pro- 
portional to cos30p, are often referred to as trigonal 
warping. Recall now this expansion is valid if the con- 
dition IpI ^ h/a is fulfilled. Rewriting the expression 
(z;jr/6i)|p| = IpI a/ (4/1), higher order terms in can 
thus be viewed as a perturbation of the original eigen- 
value = ±Vf\p\ in the small parameter (Aw^^lpl), 



where A = a/6t will be used below to identify the order 
in the perturbation. In the semiclassical limit (fi, — > 0), 
only the modification of the action needs to be taken into 
account since this latter is multiplied by the large param- 
eter 1/h. Our aim is therefore to compute the (first and 
second order here) corrections to the action in an expan- 
sion in A 



(51) 



In the presence of a constant magnetic field, the clas- 
sical equations of motion derived from the first order 
approximation = ±Ujr|n| are integrable, and this 
property is not modified by the addition of terms in 
depending only on |n|. This can be easily shown by per- 
forming a canonical transformation to the guiding cen- 
ter coordinates. For sake of completeness, this canonical 
change of variables is detailed in appendix |BJ The new 
coordinates read 

R- = (^n,y,a;o) 

P = (n^eByo) 

with ro the center of the cyclotron orbit, so that |n| = 
y^Pf+lTBXy and tan(/.n = eBX/P^. We thus have 

H+ = -■it' + p- pL2\p^ - MiA^^ 

with Vp[Px + ieBX) — pe"^, ^2 — (cosScj) — Qa^), and 
^1 = 1(3 -f cos^ 3(f> + 6a Y cos 30). In this new system of 
coordinates, the action is easily calculated as 

with the constraint E — Therefore, to order A^, and 
with E' = E + 3t' 

p2 ^ ^'2 + 2^iXE'^ + (5/ii^ + 2pL2)X^E"^ 

which gives for the action 

1 / i' 

S = — 2ttE'^ - 2ATTa-XE'^ 

2vleB \ t 

+ 127r(l + 30(y)2)A2s'4 

The third order terms had to be taken into account in 
the low-energy expansion, since their contribution in the 
second order correction of the action is of the same mag- 
nitude as that of second order terms. The third order 
term in the next-to-nearest neighbor contribution how- 
ever cancels out in the calculation of 5* and thus a second 
order expansion in h'{p) would have been sufficient. In- 
troducing this shift in the action in the Landau-levels 
calculation of section IIV Al finally gives 



= £;„ ( 1 ± Qa-XEn - iX^El 



(52) 
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(Ib = •\/fi./(ei3) is the magnetic length). This result is 
identical to the one obtained purely quantum mechan- 
ically in [23] ■ As discussed in this paper, the resulting 
effect is however too small to interpret shift in Landau 
levels observed experimentally by Plochocka et al. [13] ■ 



D. Effect of a mass term 

To end this section, let us consider the effect of a con- 
stant mass term rnQV^Oz in the graphene Hamiltonian, 
so that 



= ±Jmlv% + vlW- 



(53) 



Interestingly, a constant mass term does not modify the 
time derivative M{t) of the semiclassical Berry-like phase 
since (see Eq. ([TS])) it depends only on the gradient of 
m(r). Furthermore, as shown by a direct calculation, 
the energy dependence of the Landau frequency is not 
affected either by the mass term. Therefore 



T 

Mm) 



2-K _2'K E 
^ vleB 



OJ 



eB_ 

'2E 



const. 



and the semiclassical phase 



Mj{t)dt = jan 



(54) 
(55) 

(56) 



is the same as without the mass term. 

The mo dependence of the Landau level position is 
therefore entirely due to the uiq variation of the action 



Sj ^ jir 



eBvl 



1 



(57) 



which, following the same steps as in section HVAl gives 
p{E) = {A/{2^ll)) Y.t=, S{E ± En), with 



En — 



^El{Q) + mlv% 

1 (moWp)^ 



En{Q) 1 



An ehB 



(58) 
(59) 



(i?„(0) is the value of E„ at mo = given by Eq. (gS]))- 
One recovers semiclassically in this way the result origi- 
nally derived by Haldane |25| . 



V. SEMICLASSICAL VERSUS ADIABATIC 
BERRY PHASE 

We would like to finish this paper with some general 
discussion concerning the semiclassical phase 

U = - i MMt),r{t))dt (60) 



Af,c(p(<),r(t)) = Im 



V 



±t 



dH fdV 



9p 



r± 



dr 



(61) 



(see Eq. (O) computed on a periodic orbit (p(i),r(t)) 
(of period T). 

That, for a clean graphene monolayer without a mass 
term, ^sc = T^r (as expressed by Eq. ([38|) . with j = 1) is 
usually said to be expected since the corresponding con- 
figuration is exactly the one discussed in detail by Berry 
in his 1984 paper 9]: the path of integration corresponds 
to encircling once the Dirac point, where the and 
H~ manifolds intersect. This argument however relies 
on an exact intersection between the two manifolds, and 
should a priori not apply when a mass term mo intro- 
duces a gap. From this perspective, one does not expect 
the Berry phase to be equal to ±7r when mo ^ 0, and 
Eq. ([56]) may come asji surprise. (Note though this was 
already observed in [2(|.) 

The resolution of this apparent paradox is that, as dis- 
cussed in [l^ [l^ , the semiclassical phase ^sc defined by 
Eq. (PT|) and the adiabatic phase introduced by Berry are 
closely related, but eventually different, quantities. Both 
of them are induced by the adiabatic variation of the 
eigenstates and V~ along the trajectory. However, 
the point of view taken in the semiclassical approach is 
that both the internal space (associated here with the 
sub-lattices {A, B)) and the external space (position r) 
are coupled dynamical variables. Treating the coupling 
between these variables in the semiclassical approxima- 
tion (which indeed implicitly assumes that the "external" 
variable is slow and the internal variable fast) leads to the 
semiclassical expression (pT|) . 

Theproblem Berry was considering in his seminal ar- 
ticle Q is however different: in that case, only the in- 
ternal degree of freedom is considered a dynamical vari- 
able, and the external degrees of freedom are actually 
a space of parameters assumed to be entirely controlled 
by the experimentalist. One may in that case of course 
choose this path as the classical trajectory {r{t)) (with 
H{r) = iJ(p(r),r)) determined by the dynamics in the 
semiclassical approach. In that case however the corre- 
sponding phase is given by [§] 



bad 



A'U{r{t))dt 



Mad(r(t)) = ^V^^^■r 
or 



(62) 
(63) 



(the normalization of ensures that Mad is real). 

Let us assume, for this discussion, that we are inter- 
ested in the evolution of the eigenstate associated 
with the positive eigenvalue . Furthermore, let us 
switch to the bra/ket notation for the eigenvector and 
write = |±), V^^ = (±|. First order perturbation 
theory implies r = dH'^/dp = {-{]{dH/dp)\+), and there- 
fore Eq. ([55]) can be rewritten as 

fill fi 
Mad(r(t))=*(^^K-)-(^^>- (64) 

On the other hand, inserting the identity I2 ~ |+)(+| + 
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|-)(-| in Eq. m, 

(65) 

Thus, the adiabatic and semiclassical phases actuaUy dif- 
fer from the quantity 



•Cad Csc 



Im 



dt . 



(66) 



From this expression, we can see that in the absence of 
a mass term, but for an arbitrary electrostatic potential 
U{r), the semiclassical and Berry phases are identical. 
Indeed, for m(r) = 0, the expressions Eqs. (|A2|) - (|A3p for 
the eigenvectors of H{p,r) take the simple form 



1+^ V2 («e'"'^ 
with (j) the phase of II^ + iHy . As a consequence 



(67) 
(68) 



(M ^ l-> • H ^) = Y (- '^^-'^ + cos cj^dycj^m 

, ,dH . ^ , , d , .aVf. , , „ , . , „ 

^ I — (cos0d:E(/) + sm(/)dy0) 



la dcp 
Y~dt ' 



(70) 



The right hand side of Eq. (|69p is purely real, implying 
that, in the simple case m = considered here, ^ad ~ 
^sc = 0. Eq. ([70]) then expresses that, independently 
of the nature of the electrostatic potential U{r), the - 
here identical - Berry phase and semiclassical phase are 
just given by plus or minus (depending on a) half the 
angle of rotation of the velocity vector. In particular, as 
demonstrated by Berry from geometric arguments Q , we 
see here from a direct calculation that for a periodic orbit, 
Cad = Csc = — ajTT, with j the number of windings of the 
trajectory. This makes particularly simple the inclusion 
of the semiclassical phase in the Gutzwiller trace formula 
Eq. (|35p when m = 0. 

Similarly, for the bilayer Hamiltonian Eq. (|30p with 
m(r) = 0, we have 



1+) 



1 



(71) 





I-) = 


1 /e-'20\ 


(72) 


with (f> the phase of Il^r - 


- illy , and 






H^> - 


J-4(- sin (j)dx(l) 
m* 


+ cos 09.y(/)j(73) 






iJ — L(cos 4>9x4' ~ 
m* 


h sin (j)dy4>) 






d(j) 

i—r ■ 
dt 


(74) 



Again, the Berry phase and semiclassical phase are iden- 
tical if m(r) = (as Eq. ([73|) is purely real), and both 
phases are given by the angle of rotation of the velocity 
vector. 

For both bilayer and monolayer graphene, it has to be 
born in mind however that in the generic case m(r) 7^ 0, 
the semiclassical phase ^sc should in general differ from 
the Berry phase ^ad- Furthermore, we do not have a gen- 
eral argument constraining any of the two phases to be 
directly related to the winding of the velocity vector (be- 
yond the case where both the mass and the electrostatic 
potential are constant). 



VI. CONCLUSION 

To conclude, we have derived an expression for the 
semiclassical Green's function in graphene and discussed 
in particular the semiclassical phase associated with the 
internal pseudo-spin structure. If no mass term is in- 
cluded in the graphene Hamiltonian, this semiclassical 
phase is identical to the corresponding (adiabatic) Berry 
phase. In that case both phases are, up to a sign, given 
by half the angle of rotation of the velocity vector. For 
a bilayer of graphene, the same result holds but with a 
phase which is twice as large. 

When a mass term is introduced however, the semi- 
classical and Berry phases in general differ. In particu- 
lar, for a clean graphene sheet in a constant magnetic 
field, we have shown that the semiclassical phase re- 
mains unmodified upon the inclusion of a constant mass 
term m(r) = jtiq, while the corresponding Berry phase 
Cad = [mo/{E — Uq) — IJajTr shows some dependence on 
toq. We have shown furthermore that in this case, what 
is relevant to the calculation of the Landau levels is the 
semiclassical, rather than the Berry, phase. Other appli- 
cations of our semiclassical formalism were also discussed, 
including the effect of higher order terms of the graphene 
Hamiltonian - e.g. trigonal warping - on the position of 
the Landau levels. 

The semiclassical approximation to the graphene 
Green's function should prove a useful tool when con- 
sidering confined electron systems in graphene, such as 
graphene nanoribbons, or more complicated geometries. 

We have benefited from helpful discussions with E. Bo- 
gomolny, J.-N. Fuchs, M.-O. Goerbig, G. Montambaux 
and F. Piechon, and thank as well all active participants 
of the weekly graphene "journal-club" held in LPS, Or- 
say. 



APPENDIX A: IMAGINARY PART OF THE 
OPERATOR (02) 

In this appendix, we give the details of the computa- 
tion of the imaginary part Im(n(2)) = (see eq. ((T7)) ) 
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of the operator 



□ 



(2) 



±t 



Here 



^+(p,r) 



Tx—^ ^ 

ox ay 



m{r)vl +e(p,r) 



v/2(e(p,r))(TO(r)w2 +e(p,r)) 

^^(0:11:^ - iUy) 
-{m{r)vl + e(p, r)) 

^2(e(p,r))(m(r)w2 +e(p,r)) 



(Al) 



(A2) 



(A3) 



are the normaUzed eigenvectors of the classical Hamilto- 
nian (see Eq. (O), 



e(p,r) = fl"+(p,r) - U{r) = ^/m{r)^vi + vlW^ 



and, with respect to spatial derivation, it is understood 
that F± = y±[(9S'±/9r),r]. 

We perform here the calculation for M+, the one for 
M~ being essentially identical. We have 



aim [V+^a^id^V+)] 

and in the same way 

Im [V+^ayidyV+)] 



2e{mVp 


+ e) 






2e{mVp 


+ e) 


1 




2e{mVp - 




avp 




2e{mVp - 





Im 



/ TT -TT \\ I dximvl. + e) 



[{mvl + e)dx:Ily - Ilydx:{mvl + e)] , 



[mv 

,2 I 



(A4) 



,+£,Vp(allx liiy)) ay y^^dyiaUx + iUy) 

e){--dyUx)+Il,dyimvl+e)] , (A5) 
I 



so that 

Im 0(2) = — ^ — X n + "-^-^ '- 

^ ' 2e \or mvi + e 



with the unit vector in the direction perpendicular to 
the graphene plane. 
Using finally that 



dr 



B 



{dr X drS)^_ = d,dyS-dydxS = 

and that the Hamilton-Jacobi equation E 
H+{^,r) = implies 



d ( .dS+ . 



dU 
dr 



then gives Eq. (fTSi 



Then, introducing the guiding center Tq = {xo,yo) coor- 
dinates 



yo 



1 

eB 



Py 



y 



2 eB' 

we define the point transformation R(r') as 



This transformation is obtained from the generating 
function F(r',P) = P.R(r'), and therefore the new mo- 
mentum is given by 



APPENDIX B: GUIDING CENTER 
COORDINATES 

We sketch here the construction of the new canonical 
variables (R — (AT, F),P — {Px,Py)) introduced in sec- 
tion llV C] We start first by performing the simple canon- 
ical transformation r ^ r' = (a;,py),p ^ p' = {px, ~y)- 



which is easily inverted into 



,s , eB 
P(P )^{Px- —y,Px 



eB 



y) = {\lx,eBya) 



The unperturbed Hamiltonian is then given as 

+ (p, r) = |n| = «^ ^fifTWxf. 
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